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DERIVATION OF THE LINEAR LANDAU EQUATION AND LINEAR 
BOLTZMANN EQUATION FROM THE LORENTZ MODEL WITH MAGNETIC 

FIELD 


M. MARCOZZI AND A. NOTA 


Abstract. We consider a test particle moving in a random distribution of obstacles in the plane, under 
the action of a uniform magnetic field, orthogonal to the plane. We show that, in a weak coupling limit, 
the particle distribution behaves according to the linear Landau equation with a magnetic transport 
term. Moreover, we show that, in a low density regime, when each obstacle generates an inverse power 
law potential, the particle distribution behaves according to the linear Boltzmann equation with a 
magnetic transport term. We provide an explicit control of the error in the kinetic limit by estimating 
the contributions of the configurations which prevent the Markovianity. We compare these results with 
those ones obtained for a system of hard disks in IBMHHj . which show instead that the memory effects 
are not negligible in the Boltzmann-Grad limit. 


Keywords: Lorentz gas; magnetic field; linear Boltzmann equation; linear Landau equation; low density 
limit; weak coupling limit. 
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1. Introduction 


Consider a point particle of mass to = 1 in d = 2, 3 moving in a random distribution of fixed 
scatterers, whose centers are denoted by (ci,..., C/v)- 

We assume that the scatterers are distributed according to a Poisson distribution of parameter /i > 0. 
The equations of motion are 


( 1 . 1 ) 



here (a;, v ) denote position and velocity of the test particle, t the time and A = for any time dependent 
variable A. 


To outline a kinetic behavior it is usually introduced a scaling of the space-time variables and the 
density of the scatterer distribution. For this model, it is more physically intuitive to transfer the scaling 
to the background medium. More precisely, let £ > 0 be a parameter indicating the ratio between 


Date: February 8, 2016. 


1 





2 


M. MARCOZZI AND A. NOTA 


the macroscopic and microscopic variables, we keep time and space fixed and rescale the range of the 
interaction and the density of the scatterers, i.e. 


( 1 . 2 ) 


= £ a (j){ f) 

[J, e = ^£ _ ( d_1+2 “) 


where d = 2,3 is the dimension of the physical space and a £ [0, is a suitable parameter. This means 
that the probability of finding N obstacles in a bounded measurable set A C R d is given by 


(1.3) 


,N 


t (dc N ) = e ^ dci ■ ■ ■ dcN 


where cjv = ci,..., Cn and |A| = meas(A). Consequently, the equation of motion (1.11 becomes 


(1.4) 


V = —£ 


E,v</>(^) 


Now let T* N ( x , v ) be the Hamiltonian flow solution to equation ( |1.4| ) with initial datum (x, v ) in a given 
sample of obstacles (skipping the £ dependence for notational simplicity) and, for a given probability 
distribution fo = fo(x,v ), consider the quantity 


(1.5) 


fe(x,v,t) =E e [fo(T c *(x,v))} 


where E e is the expectation with respect to the measure P e given by (1.3). 


In the limit £ -A 0 we expect that the probability distribution (1.51 will solve a linear kinetic equation 


depending on the value of a. If a = 0 the limit corresponding to such a scaling is called low-density (or 
Boltzmann-Grad) limit. In this case f e converges to the solution of a linear Boltzmann equation. See 
0 S) BBS. ITT' . On the other hand, if a = ! 2 the corresponding limit, called weak-coupling limit, yields 
the linear Landau equation, as proven in [KPl iDGLl IK]. The intermediate scaling, namely a £ (0, |), 
although refers to a low-density regime, leads to the linear Landau equation again, see [dHE]. 

We want to remark that in iDPlIDRj the authors exploit the original constructive idea due to Gallavotti 
(see m for the Boltzmann-Grad limit. This method is based on a suitable change of variables which 
can be implemented outside a set E of pathological events which prevent the Markov property of the 
limit (such as the set of configurations yielding recollisions, i.e. when the test particle recollides with a 
given obstacle after having suffered collisions with other different obstacles). The probability P e (-E) is 
vanishing as £ tends to 0. The main difference is that in m the range of the potential is infinite in 
the limit, therefore the test particle interacts with infinitely many obstacles. As for the case of the long 
range potential considered in m also in |DR| there is a lack of the semi explicit form of the solution of 
the limit equation. This requires explicit estimates for the set of bad configurations of obstacles. For a 
short range potential, like in the case of the hard-sphere potential considered in 0, a simple dimensional 
argument is sufficient. For an explicit control of the error in the kinetic limit for the liard-sphere potential 
see for instance [ BNPP1 . Moreover, in [DRllK] it was proven that even if a > 0, but sufficiently small, the 
recollisions are still negligible. Incidentally we note that, if a is close to 1/2, this is not true anymore and 
it would be interesting to derive the Landau equation in this regime, by means of an explicit constructive 
approach. 

Furthermore it has been observed that the presence of a given external field, in the two dimensional 
Lorentz model, strongly affects the derivation of the linear Boltzmann equation in the Boltzmann-Grad 
limit. Bobylev et al, in [BMHHj and later in 1BM11H ll 1BHPH] (see also KSl for further readings), showed 
that the set of pathological configurations is no longer negligible when the test particle moves in a plane 
with a Poisson distribution of hard disks and a uniform and constant magnetic field perpendicular to the 
plane. See Figure [l] for a pictorial representation of the light particle’s motion. 

The following simple computation turns out to give a good heuristic argument explaining these results: 
consider the probability of performing an entire Larmor circle without hitting any obstacle, Rl being 
the Larmor radius. From equation (1.3) one easily gets 


Pfl; ~ e ~^ eArea (^e) ^ g-2 

where £/ e (Rl) is the annulus of radius Rl and width e. Hence, P# t is not vanishing in the limit £ —> 0 
and the Markovianity of the limit system can not be attained. In fact, in | BMHH1« IBHPH] , a kinetic 
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Figure 1. Typical paths of the charged test particle when the obstacles are hard disks: 
due to the magnetic field, it performs arcs of circle between two consecutive collisions. 


equation with memory is derived, i.e. a generalized Boltzmann equation, taking into account those effects: 


( 1 . 6 ) 


D 

Dt 


[t/T L 


f G (x,V,t) =H e £ ^ 


a -vkT L 


k —0 


IS i 


dn (v ■ n) 


[\{v ■ n)b n + x(~v ■ n)]/ ( x, S 0 z v,t- kT L ), 


where f(x,v,t) is the probability density of finding the moving particle at time t at position x with 
velocity v and 


(1.7) 


f G (x,v,t) 


f{x,v,t ) if 0 < t < Tl 

(1 - e~ vTL )f{x,v,t) if t>T L . 


Here v = 2\v\ p e e is the collision frequency and Tl = 2ir/£l is the cyclotron period where f 1 = qB/m is 
the frequency, being q the charge and m the mass. Furthermore, note that 

— = (d t + v • V x + (v x B) ■ V„) 

is the generator of the free cyclotron motion with frequency f1 and [; t/Ti,\ the number of cyclotron periods 
Tl completed before time t. The angular integration over the unit vector n in (1.6) is over the entire unit 
sphere Si centered at the origin. In the gain term the operator b n is defined by 


b n <p(v) = <f>(v — 2(v ■ n)n) 

where (f>(v) is an arbitrary function of v. The precollisional velocity v' = v — 2{v ■ n)n becomes v after 
the elastic collision with the hard disk. Note that v' ■ n < 0. In the loss term, the precollisional velocity 
v is also from the hemisphere v ■ n < 0. Finally, the shift operator S^ k , when acting on v, rotates the 
velocity through the angle —k6 1 where 6 is the scattering angle (from v' to v). 

For further readings in this direction we refer to [DR111DR2I , where the authors consider a stochastic 
Lorentz model with a smooth external force field F(x, t) and with absorbing obstacles, i.e. the interaction 
between the obstacles and the test particle is such that the test particle disappears whenever it enters 
an obstacle. It is proved that the kinetic equation associated to this model in the Boltzmann-Grad limit 
is non-Markovian and that the Markovianity can be recovered by introducing an additional stochasticity 
in the velocity distribution of the obstacles. 
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In this paper we consider the case of a random distribution of scatterers in R 2 where each obstacle 
generates a smooth positive and short-range potential <f>, with a > 0 and sufficiently small. We show 
that, in this case, the solution of the microscopic dynamics converges, in the intermediate limit (when 
a £ (0,1/8)), to the solution of the linear Landau equation with an additional transport term due to the 
magnetic field. From the heuristic point of view, this result is suggested by the observation that in this 
case the probability F Rl of performing an entire Larmor circle without hitting any obstacle is given by 

^ e ~^ e 2neR L ^ e -2-KR L fie~ 2o ‘ 

which vanishes as e —> 0. This computation shows that one family of the pathological events preventing 
the Markovianity is negligible in this setting. We stress that this rough argument is not sufficient to 
conclude that we can recover the Markovianity in the limit. Indeed, to prove this, we need to show that 
all the other bad configurations of obstacles defining the set E are negligible in the limit, as we will see 
in Section 14711 

Furthermore, we observe that even if we consider a long range inverse power law interaction potential, 
truncated at distance e 7 ” 1 with 7 £ ( 0 , 1 ) suitably large, in the low density regime a = 0 , we can prove 
that the memory is lost in the limit. More precisely, we prove that the microscopic solution converges 
to the solution of the uncutoffed linear Boltzmann equation with a magnetic transport term. With the 
same purpose of the rough argument presented above, we observe that the probability P r l of performing 
a complete Larmor circle without hitting any obstacle is approximatively given by 

F Rl ~ e ~2irR L ^s' y ~ 1 

which vanishes as e —» 0 when 7 < 1 . Also in this case this represents only one example of bad 
configuration of scatterers. It is essential to prove that the contribution of the whole set of pathological 
events is negligible in the limit, as we will show in Section[5] Moreover, from the technical point of view, 
we observe that the parameter 7 has to be chosen close to 1 as dictated by the explicit control of the 
memory effects. 

Thus, as we pointed out with the heuristic motivations above, the non-Markovian behaviour of the 
limit process, discussed in |BMHHj , disappears as soon as we slightly modify the microscopic model given 
by the two dimensional Lorentz Gas. 

The purpose of this paper is to provide a rigorous validation of the linear Landau equation and the 
linear Boltzmann equation respectively with magnetic field by using the constructive strategy due to 
Gallavotti. We remark that, as in [DPI DRl IBNP1 iBNPPj . we need explicit estimates of the error in 
the kinetic limit and this is the crucial part. Moreover, as a future target, it could be interesting to 
understand if a rigorous derivation of the generalized Boltzmann equation proposed in IBMHHj can be 
achieved by using the same constructive techniques. 

The plan of the paper is the following: in the next Section we establish the model and formulate the 
results; in Section [3] we present the strategy of the proofs, whereas Sections [4] and [5] are dedicated to 
the nontrivial analysis and explicit estimates of the sets of bad configurations producing memory effects, 
which is the technical core of this paper. 


2. The Model and main results 


2.1. The Lorentz Model with short range interactions. We consider the system (1.4) in the plane 
(d = 2) under the action of a uniform, constant, magnetic field orthogonal to the plane. The equations 
of motion are 


x = v 

v = Bv - 1 - V(ft( |x ~ Cil ) , 

where B is the magnitude of the magnetic field and u _L = ( 1 ) 2 ,— 17 ). We assume that the potential 
<f> : R + -A R + is smooth and of range 1 i.e. 4>{r) = 0 if r > 1. Therefore the particle is influenced by the 
scatterer a if \x — Cj| < e. 

Starting from the initial position x with initial velocity v, the particle moves under the action of the 
Lorentz force Bv 1 ". Suppose that the particle has unitary mass and unitary charge, namely m,q = 1, hence 
between two consecutive scatterers, the particle moves with constant angular velocity f l = qB/m = B 
and performs an arc of circle of radius Rl = \v\/B. Rl is the Larmor radius, i.e. the radius of the 















5 


cyclotron orbit whose center is situated at the point 


x c = x + 


&{%) -v 

n ’ 


where the tensor &(g>) denotes the rotation of angle <p. Without loss of generality we assume from now 
on that |u| = 1, therefore Rl = 1/B. Moreover, we will denote by Si the kinetic energy sphere with 
unitary radius. 

The precise assumptions on the potential are the following: 

Al) 4> £ C 2 ([0, oo)); 

A2) > 0, < 0 in (0,1); 

A3) supp (j> C [0,1]. 

On /o we assume that 

A4) / 0 £ C 0 (R 2 xl 2 ) is a continuous, compactly supported initial probability density. Suppose also 
that \D k x f 0 \ < C , where D x is any partial derivative with respect to x and k = 1,2. 

Moreover, we assume that 

Ah) The scatterers are distributed according to a Poisson distribution (1.3) of intensity g e = pe~ 5 
with 5 = 1 + 2a, a £ (0, |). 

Next we define the Hamiltonian flow T* (x, v ) associated to the initial datum ( x , v), solution of (2.1) for 
a given configuration cjy of scatterers, and we set 

(2-2) f E (x,v,t) = E,[f 0 (T-*(x,v))} 

where E e denotes the expectation with respect to the Poisson distribution. 

The first result of the present paper is summarized in the following theorem. 


Theorem 2.1. Let f e be defined in (2.2). Under assumption Al) — A5), for all t £ [0,T], T > 0, 

!im =g{-,t) 

£—>•0 

where g is the unique solution to the Landau equation with magnetic field 

j (d t +v ■ V x + Bv 1 - ■ V v )g(x, v, t) = £, A Sl g(x,v,t) 

1 g(x,v, 0 ) = fo(x,v) 

where A g 1 is the Laplace-Beltrami operator on the circle S\ and £ > 0. The convergence is in L 2 ( 


(2.3) 


! x50. 


The constant ^ is the diffusion coefficient and its explicit expression will be given below in (3.9) and 
in remark f3.41 


2.2. The Lorentz Model with long range interactions. We consider now the case in which each 
obstacle generates a potential of the form 


$e{\x~c\) = —-) 

£ 


where the unrescaled potential il> e is an inverse power law potential truncated at large distances. More 
precisely we assume the following: 

Id < e 7_1 


B 1) 


ip s (x) = 


= J 


_-s( 7 -l) 


> e 7 " 1 


with 7 £ ( 0 , 1 ) and s > 2 . 


We point out that it could be challenging to consider directly the untruncated long range potential 
^>(|:r|) = |a:| _s . In fact, this problem presents deep additional difficulties as noted in Remark 2.3 in 
and new ideas and techniques are necessary. 

Moreover, we assume that 

B 2) The scatterers are distributed according to a Poisson law (1.3) of intensity g e = £“V, g > 0. 
The equation of motion in macroscopic variables reads 


x = v 
v 


l) = Bv ± -e 


(2.4) 

with given in Assumption BA). 

Let T 7 v (x.v) be the Hamiltonian flow solution to equation (2.4) with initial datum (x,v) in a given 
sample of obstacles. Let be /o = fo(x,v) be the initial probability distribution. On /o we assume 
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B3) f 0 £ L 1 D fF 1,00 (K 2 x R 2 ), /o > 0, J f 0 dxdv = 1. 

We consider the quantity 

(2-5) f E (x,v,t) = E e[fo(.T~*(x,v))} 


where E e is the expectation with respect to the measure 
paper is summarized in the following theorem. 


given by (1.3). The second result of this 


Theorem 2.2. Let f e be defined in (2.5). Under assumption B 1) — B 3) with 7 £ (6/7,1), for all 
t £ [0, T], T > 0, 

lim = f(-,t) 

E ->0 

where f is the unique solution to the linear Boltzmann equation with magnetic field 

j (d t + v - S7 X + Bv 1 - ■ V„)/(i, x, v) = Lf(t, x, v ) 

1 f(x,v, 0) = fo(x, v), 


( 2 . 6 ) 

with 


Lf(v) = H [* T(0) {m(e)v) - /(«)} d<9 , 

J —TV 


T(0) is the differential cross section associated to the long range potential ip(\x\) = \x\ s and the operator 
rotates the velocity v by the angle 0. The convergence is in 3T(M. 2 x Si). 

3. Proofs 


3.1. Proof of Theorem 2.1, Following |DP . DRl IBNPl we split the original problem into two parts. 
The first one concerns the asymptotic equivalence between f e defined in (1.5) and h e , solution of the 
following Boltzmann equation 

j (d t + v ■ X7 x + Br 1 • W v )h e (x,v,t) = L e h e (x,v,t) 

1 h s (x,v,0) = fo(x,v) 


(3.1) 

where 

(3.2) 


L E h e (v) = He J dp[h E (v') - h s (v)]. 


Here v' = v — 2(w • v)u> is the outgoing velocity after a scattering with incoming velocity v and impact 
parameter p £ [—e,e] generated by the potential e a </(j). Moreover, w = u(p) is the versor bisecting the 
angle between the incoming and outgoing velocity and 9 e is the scattering angle. The precise result is 
stated in the following proposition. 

Proposition 3.1. Under assumption HI) — H5), for any T > 0, 

(3-3) lim ||/ e - /i £ ||Loo([ 0 T ] ;i i( R 2 X 5 1 )) = 0 


where h e solves (3.1). 

The proof of the above Proposition is postponed to Section [4} 

The second step concerns the grazing collision limit. Note that the presence of the magnetic field does 
not affect the last step. More precisely we have the following 

Lemma 3.2. The deflection angle 9 e {p) of a particle colliding with impact parameter p with a scatterer 
generating a radial potential e a f> under the action of the Lorentz force Bw L satisfies 

(3.4) \6 e {p)\<Cs a . 


Proof. As established in [jDR. (Section 3), the estimate (3.4) holds when the test particle scatters with 
no external field. Hence, we just need to compare the dynamics of the test particle in presence of the 
constant magnetic field with the free dynamics. Let (x(t),v(t)) be the solution of the following 


(3.5) 


V = — £ 


a -1 


Vd-(^) 


Let r be the collision time for the dynamics described by GED- The key observation is that the presence 
of the magnetic field does not modify the estimate for the collision time related to the dynamics in (3.51. 
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Indeed also in our case r < Cs, C > 0, as in [DEI, see Appendix [A] for the detailed computations. 
Therefore we get 

\v(t)-v(t) | = 


< £ 


< £° 


a-1 


< e' 


a—2 / 


< £ 


a—2 / 


ds (F(x(s)/e) — F(x(s)/e)) + / dsv ± B 

Jo 

ds | F(x(s)/e) — F(x(s)/e) \ + C\£ 

JO 

2 C 2 f ds |x(s) — x(s)| + C\£ 

Jo 

2 C 2 f ds f dt\v{t) — v_{t)\ + C\£ 

Jo Jo 

2 C 2 f ds f dt\v(t) — y(t)\ + C\£, 

Jo Jo 

where F(x) := — fS/<f>)(x). By using Grdnwall’s inequality we obtain 

(3.6) |u(r) — v(t)\ < Ci £e Cse T <C'i£e c ' 3E 

for a > 0 and £ sufficiently small. Hence, the velocities v and v are asymptotically equivalent up to an 
error term of order £. We now define v' and yj to be the outgoing velocities with and without magnetic 
field respectively, v the incoming velocity. By using (3.6) we have 

(3.7) a/ 2(1 — cos 9 e ) = \v' — v\ < \v' — y'\ + |u 7 — v\ < Ce + \J 2(1 — cos 6 e ) 

where 9 e is the scattering angle without magnetic field. From iDR i we know that 9 e < C’E a , so from 
(3.7) we obtain 

I sin ^ | < Ce + | sin ^ | < Ce + C'E a < C"s a . 

Since 9 e is continuous as a function of the impact parameter p , it results 9 e < Ce 01 . For further details 
see Proposition |B.l | in Appendix [B] □ 

The following proposition shows the asymptotic equivalence between the solution of Landau equation 
and the solution of the previous Boltzmann equation h e . 

Proposition 3.3. Under the assumptions Al) — A5), h e —>• g in L°°([0, T]; L 2 (R 2 x Si)) where g is the 
unique solution to the Landau equation with magnetic field 

(d t + v ■ S7 X + B V 1 ■ S7 v )g(x, v,t) = £ A Sl g(x, v, t ) 
g(x,v,0) = f 0 (x,v) 


(3.8) 
where 

(3.9) J ^9 2 e (p) dp 

is the diffusion coefficient. 

Remark 3.4. As shown in Appendix |7j[ 9 e = 9 e + 0(e), where 9 e is the scattering angle without any 
magnetic field, i.e. the scattering angle studied in m- This implies that the explicit expression for the 
diffusion coefficient obtained in EE! still holds in our case: 


(3.10) 


^ 2 




where the integrand is an even function of the impact parameter p. 


Remark 3.5. The linear Landau equation (3.8) propagates the regularity of the derivatives with respect to 
the x variable thanks to the transport operator. Moreover, the presence of the collision operator A? := 
lets the solution gain regularity with respect to the transverse component of the velocity. Indeed, under 
the assumption A4) on ff, the solution g : R 2 x Si —> R + satisfies the bounds 

(3.11) \D k x g\<C, \D*g{x,v)\<C Wk<2,h>0, 

Vt £ (0, T], where C = C(fo,T) and D v is the derivative with respect to the transverse component of the 
velocity. In particular, the solutions of (3.8) we are considering are classical. 
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Proof. By using the invariance of the scattering angle with respect to the space scale, we rewrite the 


collision operator in the right hand side of (|3.l|) as 
(3.12) 


L e h E {v) = p £ e j dp[h e (v') - h £ (v)\. 


We look at the evolution of h £ — g, being g the solution of (3.81, namely 
(3.13) 


(d t +v-X7 x + By 1 ■ V„) (h e - g) = (h £ h £ 


where Jz? := £ A^. 

_Note that g £ L 2 (R 2 x Si) because fo £ L 2 (R 2 x K 2 ) and h £ £ L 2 (K 2 x Si). Indeed, from Proposition 
we know that h e £ L X (K 2 x Si) but the hypothesis on the initial state implies that h e £ L 2 (R 2 x Si). 
We now consider the scalar product of equation (| 3. 13 ) with ( h e — g) in L 2 (K 2 x 5i) and we obtain 


3.1 


-d t \\h e -g\\\ =- {h £ -g, -L e [h e - g\) 

+ ( h £ -g , [L e 

By exploiting the positivity of —L E and the Cauchy-Schwartz inequality we get 

dt\\h e -g\\ 2 < || (L e — Jf)g\\ 2 . 

We now set 

g{v')-g(v)= (v 1 - v) ■ V\ Si g{v) 

+ ^{v 1 -v)® (v' - v)V ls V l3i g(v) 

+ ^(v' -v)® (v' -v)® (v' - v)X7 ls V lsi X7 lsi g(v) + R e , 
with R e = 0(\v — i/| 4 ). Integrating with respect to p and using symmetry arguments we obtain 

Leff = ps~ 2a \\ A s i9 J dp\v'-v\ 2 + J dpR s |. 

Observe that \v' — v\ 2 = 4 sin 2 , then by direct computation 

iTo ^~2 /_, dp W ~ " |2 = ^~2 /_, ^ [p) dp =■ t - 

Therefore, thanks to Lemma [372] , we have 

||(L E -^) ff || L2 < e 2 “||Af Si5 || i2 < £ 2 “C, 


which vanishes for e — > 0. 


□ 


dp 


Remark 3.6. We avoided introducing the cross-section r(0 E ) := of the problem because the map 
p —> 9 e (p) is not monotonic in general. 

Indeed if <j> is bounded and e sufficiently small, ^v 2 > £ a (f>(0) so that 9 = 0 for p = 0 and p = ±1. ^4s a 
consequence, r(0 £ ) is neither single valued nor bounded. 


3.2. Proof of Theorem 2.2, The general structure of the proof follows the lines of [DP] where an 
analogous result has been proven when the magnetic field is zero. 


Proposition 3.7. Let f £ be defined in (2.5). Then, for any T > 0, 
(3.14) 


ll/e — ^e ,7 ||l°°([OjTJjL 1 (R 2 XSi)) — 0 


£—>-0 


where /i Ej7 is the unique solution of the truncated linear Boltzmann equation with magnetic field 

(dt + v ■ V x + Bv 1 ■ V„)/i Ej7 (f, x, v) = Lh e ^(t,x,v) 
h £ ^(x,v, 0) = f 0 (x,v) , 


(3.15) 


Lf(v) = p / P< B >) {f(0(6)v) - f(v)} d9. 


with 
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/ D\ 

and Te i7 ; is the differential cross section associated to the unrescaled potential ip e with magnetic field. 
The proof of Proposition |3.7| is in Section [5j 

This allows to reduce the problem of the transition from the solution of the truncated linear Boltzmann 
equation to the solution of the untruncated linear Boltzmann equation to a partial differential equation 
problem. Indeed, as in [DPj . we can prove the following 


Proposition 3.8. Let h en solution of (3.15). Then, for any T > 0, 


(3.16) 


where f is the unique solution of (2.6). 


h E)1 -> / in C([0, T]\ $>') 


Proof. In Appendix IbI Proposition B.2 the cross section rl[ 7 (0) is shown to be bounded by C8 1 1 t s 


and to converge to T(0) almost everywhere as e —> 0, where T(0) is the cross section associated to the 
truly long range potential 'P(r) = r~ s without magnetic field. Therefore, the proof of Proposition 3.8 is 
exactly the same as the one of Proposition A.2 in [DP]. □ 

4. Proof of Proposition I3JJ 


In this Section we prove the asymptotic equivalence of f E , defined by (2.2), and h e , solution of the 
linear Boltzmann equation (3.1), that we recall here for the sake of clarity 

(4.1) (d t + v ■ V x + B u 1 • S7 v )h e {x, v, t) = L E h e {x, v, t), 

where 


(4.2) 


L E h E (v) = pe 


dp{h e {v') - h E (v)}. 


/-1 


This allows to reduce the problem to the analysis of a Markov process which is an easier task. Indeed, 
the series expansion defining h e (obtained perturbing around the loss term) reads as 


(4.3) 


h e (x,v,t) = e 


— 2e~ 


* fit 


Ea*. 

Q> 0 


Q 


dt Q ... dt i / dpi... dp Q /o(q t {x,v)). 


Here 7 _t (x,i>) = (£ e {—t),r] E (—t)) where r] e is an autonomous jump process and f E is an additive func¬ 
tional of r] E . Equation (4.3) is an evolution equation for the probability density associated to a particle 
performing random jumps in the velocity space at random Markov times. 


We start the proof by looking at the microscopic solution f E defined by (1.5). For (x,v) £ 
t > 0, we have 

N r 

0 -ii e \3S(x,t)\ y- 

^ N\ 

N> 0 


(4.4) 


f e (x,v,t) = e~ 


!sg(x,t) N 


dc N f 0 (T*(x,v)), 


where T* (x, v ) is the Hamiltonian flow with initial datum (x, v ) and ^(x, t ) is the disk centered in x 
with radius t. 

Given the configuration of obstacles cjv = ci ... Cn, we shall say that Ci is internal if it influences the 
motion up to the time t, i.e. 

(4.5) 

while we shall call Cj external if 

(4.6) 


inf |x e (-s) - Cj| < e, 

o <s<t 


inf |x e (-s) - Cj| > e. 

0<s<t 

Here (x e (-s), v e (-s)) = T~ n s (x,v ), s £ [0, t]. 

We can perform the integration over the external obstacles and we get 

(4.7) f E (x,v,t) = [ 

Q'- Jss(x,t)Q 

where x{E) is the characteristic function of the event E and ^(b q) is the tube 

(4.8) Sf(fa q) = {y £ @(x,t) s.t. Us e (0, t) s.t. | y - x e (-s)| < e}. 

Note that in the previous integration we are not considering possible overlappings of obstacles. This is 
legitimate because we shall see that this event is negligible as £ tends to 0. 


db Q e ^ |,Sr(bQ)l x({bQ internal})/ 0 (T bQ t (x,u)), 
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Figure 2. On the left a cyclotron orbit completed without suffering collisions is repre¬ 
sented. On the right there is a repeated collision with the same scatterer. 


Furthermore, let us restrict to the configurations such that the light particle’s trajectory does not start 
from inside an obstacle and does not end inside an obstacle: in formula 

(4.9) Xi(bg) = x{bg s.t. bi g &(x,e) and h 3§{x E (-t),e) for all* = 1,... ,Q}. 

As for the overlappings, this choice is not really restrictive because the contribution related to 1 — Xi is 
going to vanish in the limit, as we shall see. Moreover, we will now list other events that will turn out to 
be negligible as e approaches 0. 

i) Complete cyclotronic orbit. 

A first cyclotron orbit is completed without suffering any collisions and a repeated collision occurs 
with the same scatterer without any collision in the meantime. 

We set 


(4.10) Xcirc(bQ) = x ({t>Q s.t. i ) is realized}). 

For a pictorial representation of the event i) see Figure [2] We now define 


(4.11) 


f e (x,v,t) f db Q e ^ | ^ (b « ) l x ({b g internal}) 

q> o dss { x,t)Q 

- Xdrc(b Q )) Xi(b Q )fo(T-*(x,v)). 


Note that f e > f e . For t <Tl one expects that the approximation with the dynamics of the test particle 
in absence of the external field is true. The unexpected fact is that even for t > 0 (Tl ) this still holds 
because (4.10) tends to 0 as e — > 0. Hence, for a given configuration b q such that xi[l — Xcirc](bq) = 1, 
we have that the measure of the tube can be estimated by 


(4.12) 


|^(b(j)| < 2 et. 


At this point we define 

fe(x,v,t) =e~ 


Q 


(4.13) 


n i 


Q> 0 


Q- JS8{x,t)Q 


dh Q x({b Q internal}) 


Thanks to (4.12) we get 
(4.14) 


(1 - Xarc(h Q )) xi(bQ)/ 0 (T b Hx, v)). 


fe> fe> fe- 


According to a classical argument introduced in m (see also [DPI EH [BNP ) , we now want to remove 
from f e all the events that prevent the light particle’s trajectory to be the Markov process described by 

h e - 
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For any fixed initial condition (x,v) we order the obstacles according to the scattering 

sequence. Let pi and t, be the impact parameter and the backwards entrance time of the light particle 
in the protection disk around bi , namely &(bi,e). Then we perform the following change of variables 

(4.15) bi,... ,b N -)• ... ,p N ,tN 

with 

0 < tjv < fjv—l < • • • < ti < t. 

Conversely, hxed the impact parameters {p{\ and the hitting times {U} we construct the centers of the 
obstacles bi = b(pi,ti) and a trajectory j~ s (x,v) := (£ E (—s),fj E (—s)), s £ [0,t] inductively. 

Suppose that we are able to define the obstacles bi, ..., 6j_i and a trajectory j~ s (x, v) := (£ e (—s), fj E (—s )) 

up to the time s = t*_i. We then define the trajectory between times t,_i and ti as that of the evolution 

of a particle moving under the action of the Lorentz force and of the potential £ a <^>(£ _1 | ■ —i>*_i|) with 
initial datum (£ e (— fj e {—ti-i))- Then bi is defined to be the only point at distance e of £ e (— ti) and 
algebraic distance pi from the straight line which is tangent to the trajectory at the point £ e (— ii). 

However, l~ s (x,v) = (x E (—s),v e (—s)) (therefore the mapping ( 4.15| ) is one-to-one) only outside the 
following pathological situations (relative to the backward trajectory). 

ii) Overlapping. 

If bi and bj are both internal then £8(bi, e) D t${bj,e) ^ 0. 

iii) Recollisions. 

There exists bi such that for s £ (tj+i,tj), j > i, £ E (—s) £ 33(bi,e). 

iv) Interferences. 

There exists bi such that £ e (—s) £ 38{bj,e) for s £ (U+i,ti), j > i. 

We simply skip such events by setting 

Xov = x({t>Q s.t. ii) is realized}), 

Xrec = X'({ b Q S.t. iii) is realized}), 

Xint = x({ b Q s.t. iv) is realized}), 


and defining 

(4.16) 

Note that 

(4.17) 


fe(x,v,t) =e 


-2e~ 


* fit 


E^ 

Q> 0 




C^Q-l 


dtc 


dpi 


dp Q 


Xl (1 ^ Xcircj (1 - Xov)0- - Xrec) (1 - Xint)Ml \x,v)). 


fe< fe< fe< fe- 


Note also that in (4.16) we have used the change of variables (4.15) for which, outside the pathological 
sets i), ii), iii), iv) l^ t (x,v) = (x E {-t),v e (-t)). 


Next we remove xi(l - Xdrc)( 1 - Xov)( 1 - Xrec)(l - Xint ) by setting 


h E (x,v,t) = e 


— 2e~ 


* fit 


(4.18) 


E 

Q> 0 


dti 


ptQ -1 


dtr 


dpi 


dpQfoil *(a;,u)). 


We observe that 

(4.19) 1 - Xl(l - Xov)( 1 - Xcir) (1 - Xrec) (1 - Xint) < (1 “ Xl) + Xov + Xc 

Then by (4.16) and ( |4.18 1 we obtain 

I h E (t) ~ f e (t )| < <pi(e,t) 

with 

riQ -1 r e 


(4.20) 


= ||/o||oo e 


—2e~ 


fit 


Em 

Q> 0 


Q 


dti 


Xr 


' Xint • 


dt Q / dpi... dp Q 


[(1-Xi) + Xov T Xcir “1“ Xrec H - Xint]- 

We can prove that tp\ is negligible in the limit. The precise statement follows. 
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Proposition 4.1. Let ipi(e,t) be defined as in (4.20). For any t £ [0, T] 

\\<Pi(e,t)\\ L i 0 

as e —> 0. 

Proof. See Section [4T] 


□ 


To complete the proof of Theorem 2.1 we still need to show the asymptotic equivalence of h E and h E . 
Notice that h e is given by (4.3) where the trajectory r )~ t {x,v) = (£ E (—t),r] E (—t)) is a jump process in 


the velocity space, i.e. the changes of velocity are instantaneous. We compare the trajectory ^y~ t (x,v) = 
(£ e (— t),fj E (—t)) with 7 -t (ir,u): being the sequence of impact times and r < Ce the collision 

time, the spacial coordinates can differ only inside the interaction disk, while the velocities can differ only 
if t e (ii,ti + t). 

In formulae we have 

(4.21) | Vs(-t) - Ve(-t) | < C 2 e a x(t -ti< Ce). 

By exploiting the regularity of the initial condition fo we get 

(4-22) ~ fo{^{-t),fj E (-t))\ < C'[Qe + e a \(t — t\ < Ce)] 

which implies 

rtQ -1 


\h e (x, v, t) — h e (x, v,t)\ < C'e 2 ^ te 2 “ ^ n® f dt x ... f dt Q 

Q >0 

x J dpi ... J dpQ [Qe + e a x(t -t\< Ce)] 


(4.23) 


< C{te l - 2a + e i+a ). 


Hence we obtained lim \\h e — /i e ||l=o([o t1xm 2 xSi) = 0- We observe that the monotonicity argument behind 
£—>•0 U ’ J 7 


this strategy, see equation (4.171, the positivity of the solution h e of the Boltzmann equation and the con¬ 


servation of mass imply that / E , h e and h £ have the same asymptotic behavior in L°° ([0, T]; L 1 (R 2 x S^)) 
when e —> 0. 


4.1. Control of the pathological sets: proof of Proposition |47l] In this section we prove Proposi¬ 
tion [4T] This makes rigorous the claim of the heuristic argument presented in the paper’s introduction. 
For any measurable function u of the backward Markov process (£ E , rj E ) we set 


[«] = e 2/ie£t ^2 P? [ dti... j dt Q 

Q >0 

J dpi...J dp Q u(f, E ,r] E ). 


Recalling (4.19) we have 

<Pi(e,t) < H/ollooEi.vKl - Xi) + Xov H - Xcir “1“ Xrec + X^ni] ■ 


We can skip the estimates of the first two contributions, i.e. E Xj0 [(l—Xi)] and ^ x ,v[Xov\, since the presence 
of the external field does not affect the classical arguments which can be found in [BNP1 IDR1 IDPj . 
However, the presence of the magnetic field and consequently the circular motion of the test particle 
strongly affects the explicit estimates of the pathological events ii), iii). Therefore, we need a detailed 
analysis for Xdr , Xrec and Xint- 

For what concerns the pathological event due to a recollision with the same scatterer (see Figure |3|, 
we observe that Xdr = 1 if there exists an entrance time tj such that \ti — ti+ 1| > Tl — r > Tj, — Ce 
for some i = 0,... Q — 1. Moreover, Xdrc = 1 also when a test particle performs an entire Larmor orbit 
without colliding with any obstacles. As explained in the introduction, the probability of this event is 
bounded from above by Cexp(—^f^£ -2a ) := c a (e). Therefore, it results 
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_ pi piQ — l p£ p£ 

^x,v[Xc ire] < C a (e) + e _2pe£t ^2^9 dt-\ ... dt Q / dp-, ... dp Q ^ \{U < U-i 

Q> 1 J~ e •* ~ £ »=1 

(4.24) < c a (e) + e _2 * Me£ ^(2p £ e) Q f dti--- f dt Q ^ x(U < U-i ~T L + Ce) 


Q>1 

where p e = /z£ -1-2q with a € (0,1/2) and t 0 = t. 
We set 


then 


Ji := 


:= f dti ■■■ f dt Q x{ti < U -1 - T L + Ce) 

Jo Jo 


= dti--- dtixiU < ti_i - T l + Ce) 


+ Q-i 


= dti 


dti -1 


(Q-i)'- 

(ti- v - T l + Ce)®-^ 
(Q-(i-l))! 


(4.25) 


Ex.wIXcirc] < c(e) + e 2t/ie£ ^ (2/z £ £) Q ^ /*. 

Q>1 


i=l 


Note that 


I-i+i — I dti • • • 


= dti--- 


dti 


dti— l 


{U ~T l + Ce )«"* 

(Q-i)! 

(ti_i - T l + C'e) c ?-( i - 1 ) - {-T l + Cs)^" 1 ) 

(Q- (*-!))! 


= I, - 


(- T l + Ce)Q~ i+1 f" 1 


(Q-i + 1)! (* — 1)! 


— Ii ~ fi—1 


T l + Ce) 
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where fi := ^ Tl (q C 1^ —||yj- We now look at 
Q Q -2 

'52ii = Qh-'52(Q-i-j)f j 

i=l j= 0 

Q-1 

(t-T L +Cef-^(Q-mQ-j) 

3=0 


< 


(Q-1)! 

1 

(Q^ij! 


{-T L + Ce)Q-iV , Y^ {Q _ l y{-T L +Ce) Q -^P 


(Q-j)! j! 


2(i - T L + Ce) q + (T l - Ce)(i - T L + Ce)^ 1 


(4.26) 


3=0 


(Q-j)! j! 


Finally we got 

^x,v[Xcirc] < C a (e) + e _2 ‘ Me£ ^ (2/x e £) Q 


Q> 1 


2(t-T L +C£)Q (T L -C'e)(t-T L + C'£) c 3- 1 

(Q-i)! (Q-i)! 

^ (t - T* + Ck)^- 1 


= c«(e) + 2e~ 2t ^ £ (2 /i E £)(i - Ti + Ce) ^ (2 Me £)« 

Q>1 

+ e~ 2tlJ ' e£ (T L - Ce)(2^e) £ (2 / z E £)«- 1 ~ Ti + 

Q> 1 


(Q-1)! 


(4.27) 


(Q-1)! 

= Ca(£) + 2e- 2{ - TL ~ Ce)lie ~' 2a (213£~ 2a )(t -T l + Ce) + e -2 _ Ce)(2n£- 2a ) 


for a > 0 and £ sufficiently small. Hence, E I t [x c irc] vanishes as £ —?> 0. 

We now consider a generalization of Xdr'■ let be Xarc the characteristic function of the event such that 
the light particle does not hit any obstacles in a time interval equal to Tl £ v , 0 < v < 1. More precisely 
Xarc = 1 if there exists an entrance time tj such that | i* — t, ;+1 | > T^e" — t > T^e v — Ce for some 
i = 0,... Q — 1. The same computations as for E x ,v[Xdrc] show that IE x ,i>[Xarc] vanishes as £ —> 0 when 
v < 2a. In other words this shows that the motion of the light particle outside the obstacles covers arcs 
of circle and corresponding angles of order at most 0(e v ). 

Next we pass to the control of the recollision event. We observe that 


(4.28) 


Xrec (l Xarc) Xrec T Xrec Xarc — (l Xarc) Xrec + Xarc 


and this implies 

®x,n[Xrec] 5s ® 2 :,r[(l — Xarc)Xrec] d~ ^x,v [Xarc] > 

but E XtV [Xarc] is vanishing in the limit £ — ► 0 as we have seen before. Therefore, we can focus on 
E Xj „[(l — Xarc)xrec\- Let ti the first time the light particle hits the Tth scatterer v~ the incoming 
velocity, vf the outgoing velocity (with respect to the backwards trajectory) and tf the exit time. 
Moreover, we fix the axis in such a way that vf is parallel to the x axis. We have 


(4.29) 


Q 

Xrec(l - Xarc) < l 1 - xirl) 

*=1 j> 1 


y’,3 

Area 


where Xrec = 1 if and only if 6j (constructed via the sequence ti, pi,... , tj, pi) is recollided in the time 
interval (tj,tj- 1). Note indeed that a recollision can occur only if the rotation angle | J2h=i+i(@ h +W»)I > 
7r where iph is the angle covered outside the obstacles in the time interval (th+i,t^), being Oh the h -th 
scattering angle. The constraint (l — Xarc ) implies that \<fh\ < C'e v . 

Hence, since \0 h + ip h \ < Cs a + C' e v < C"e v / 2 , in order to have a recollision there must be an 
intermediate velocity Vk, k = i + 1,... ,j — 1 such that 


(4.30) 


3+-vf\<CE v/2 , 


namely is almost orthogonal to Vj~ (see Figure |4j). 
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\ a-tk+i) 



Figure 4. Backward recollision. 


Then 

(4.31) 


Xrec(l - Xirc) < i 1 “ Xarc) E E E 

i — 1 j=i+2 k=i +1 


ydi.fc 
A.rec 5 


where Xrec = 1 if and only if Xr£ c = 1 an d (4.301 is fulfilled. Following [BNP] . we fix all the parameters 
Pi,. ■ ■, PQ,ti,... ,tQ but tk+ 1- The two branches of the trajectory I 2 are rigid so that, when a recollision 
occurs, the integration domain with respect to tk +1 is restricted to a time interval bounded by 

2e 

_ < 4 p 

cos Ce v ! 2 

Performing all the other integrations and summing over i,j,k we obtain 

Q Q .7-1 

(i-fa)EE E 


E„. 


i=l j=l fc=i+l 


(4.32) 


< Cse- 2 ^ 2at J2(Q~ m ~ 2)(Q - 

Q>3 ^ '* 

< C"t 3 e 1_8a , 


which tends to 0 as e goes to 0 for a. < 1/8. 

Following the strategy used in [BNPli , since a backward interference is a forward recollision, the estimate 
for the interference event can be handled by using the Liouville Theorem. 


5. Proof of Proposition 13.71 


Our aim is to prove the asymptotic equivalence of / e , defined by (2.5), and h Sjl solution of the linear 
Boltzmann equation ( 3. 15| ) , namely 

(<9 t + v ■ V x + B v 1 - ■ V v )h en (t, x , v ) = L h e ^(t, x , v ) 
h Etl (x,v,0) = f 0 {x,v) , 


which reads 

(5.1) h e ^(x, v, t) = e~ 2>J,e ' 1 1 


E^ Q 

Q> 0 


rt 2 


dt Q ... dti / dpi.../ dp Q f 0 (£ e (-t),r] e (-t)). 

















16 


M. MARCOZZI AND A. NOTA 


We recall that we can expand f E as follows: 


(5.2) 


N r 

f e (x,v,t) = e-“' 1^1 Y.% 

N>0 1 ■ JSt 


’SS{x,t) N 


dc N fo(T c * 


(x,v)), 


where T* ( x , v ) is the Hamiltonian flow with initial datum (x, v) and &(x, t ) is the disk of center x and 
radius t. We observe that the proof follows the same strategy of Proposition 3.1 (see also Section 3 in 
EE]). As before the hardest part is the estimate of the non-Markovian contribution which is summarized 
in the following proposition. 


Proposition 5.1. Let ipi(e,t) be defined as in (4.201 with the only difference that the radius of the 
obstacles is now e 1 instead of e and the collision time r is bounded by Ce 1 instead of Ce. Then for any 
t S [0,T] 

\\ipi{e,t)\\ L i -> 0 

os £ —> 0. 


Proof. Also in this case we can skip the estimates of the contributions E^ „[(1 — %i)] and E x ,v[Xov) since 
the presence of the external field does not affect the bounds in m ■ Moreover, as in proposition |4.1[ if 
we know that E x ,v[Xrec] is negligible, then the Liouville theorem guarantees that also E x ,v\Xint\ can be 
disregarded in the limit. Hence, it suffices to focus on E x ,v[Xdrc] and E x ,v[Xrec\- So we look at 

rt Pq-i ®~f 

(5.3) E X , v [xcirc\ < c 7 (e) + e“ 2 * MeE ^ ^ (2/o e e 7 ) Q dti - dt Q ^ x(U+i <U-T L + Ce 1 ) 

Jo Jo ; _n 


Q> 1 


i=0 


where c 7 (e) = C exp(—^^e 7 l ), p e = pe 1 and 7 £ (0,1). Following the same strategy as in proposition 


14.11 we obtain 


^xAXdrc] < c 7 (e) + e 2t/ieeT y^(2/z E e 7 ) Q 

Q> 1 


2 (t -T l + Ce 1 )® (T l - CgT)(t -T l + Ce 1 )®- 1 


( 0 - 1 )! 


( 0 - 1 )! 


(5.4) 


= c 7 (e) + 2e- 2(Ti - CeT) ' i£T 1 {2ne 1 - 1 )(t-T L + Ce 1 ) + e~ 2 ^- CeJ) ^ 1 {2pe 1 ~ 1 )(T L - Ce 1 ) 


which vanishes as e —> 0 and 7 £ (0,1). 

To control the recollision event we c an follow the strategy used in Section 4.1 and in EE], Proposition 
3.1. More precisely, as in Section 


4.1 


we introduce Xarc such that Xarc = 1 if there exists an entrance 
time ti such that |t$ — f i+1 | > T k /M — r > T^/M — Ce 1 for some i = 0,... Q — 1 where M is a finite 
constant and M > 1. One can easily see that E^ „[xarc ] vanishes as e —> 0 when 7 < 1. Furthermore, 

^x,v[Xrec] < E*,„[(l - xi^AXrecl+^x.viX^J] 

but E X,v[x^rc\ is vanishing in the limit £ -A 0 as we have seen before. Therefore, we can focus on 
E a| „[(l — Xarc)Xrec]■ We now distinguish the collisions as 


(5.5) 


Q F s 

(1 - X ( c^c)Xrec < (1 - Xa xj) XrL x( Sma jk <—, Vfc = *,..., j - l) 

i=l j> 1 


Q 


+ (1 - X ( arc) H Xrec x{ sma jk > — , for some k = i, ..., j - l) 
*=1 j > 1 


where Xrlc = 1 if and only if bi (constructed via the sequence ... ,ti,pf) is recollided in the time 

interval (tj,tj- 1) and a.j k (with i < k < j) is the absolute value of the sum of the angles between the 
outgoing velocity v£ from the fc-th obstacle and the recolliding velocity v~, i.e. 


j 

Otjk = Wk\ + 22 \ 6r + Tr\ 

r=k +1 

where 9 r is the deflection angle due to the r-th scatterer and ip r is the angle covered in the time interval 
(t r+ i,tf) outside the scatterers and <pj is the angle covered between the j-th obstacle before recolliding 
with bi. Here S > 0 is a suitable parameter that we will fix later. Note that, thanks to 1 — Xarc , we have 
| ip r | < 2n/M for any r. 
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As noticed in m, the constraint sino,*, < £ S /4 implies that \O r +ipr — tt| < e s for some r = i,... ,j — 1, 


thus we get 


Q r s 

-X { ™c) ^2^2xrecx{sma jk <—, Vfc = j- l)] 

i=l j> 1 


Q>0 i=l j=i+ 2 k=i 


[ dti 

r*~ 

7 

1 0 

/ 0 



tQ-i 


dt Q 


(5.6) 

/_ £ 7 

Now we note that for M big enough 


/ , £7 r e~Y 

dpi dp 2 --- dp Q x(\dk + Wk - 7r| < e 5 ). 
-e'y J— e'y J—e'l 


C £l f n 

/ dp k x(| 0 k + <fik - tt| < e 5 ) = / dd k f^(9 k )x(\0 k +Vk~A< g<5 ) 

J-e'y J —7r 

= e f dd k T^(9 k )x{\0 k + tp k -Tr\ < s s ) 

J —7r 


r 7 r ( 1 - B ’)+ e5 


= £ 


(5.7) 


< Ce 1 ^ 


d0 fc r( s 7 )(0 fc ) 


where Ti/lj (9 k ) is the differential cross section associated to the rescaled potential 7e> while F i B J(9 k ) is 
differential cross section associated to the unrescaled potential ip e . In the last line of (5.7) we used that 
(9) is uniformly bounded in e when 9 is far from 0, as shown in Appendix [b| 

Then from (5.6) one gets 

q s 

E *,» [(! - xifc } ) J2 Y XrLx( sin a jk <j,Vk = i,...,j- l)] 

*=i j >i 


(5.8) 


< 


Ce ~2tfise-' ^ < C£ 2l ~ 2+S . 


Q> 1 


For what concerns the second term in (5.5) we note that, once we fix all the variables {te}f =1 and 
{pt\® = i except tfc, by using the same geometrical argument as the one illustrated in Figure 5 of m , one 
gets that the integral over t k +i is bounded by 

2 £ 7 


(5.9) 

It follows that 


< Se 7 " 5 . 


sm ajk 


Q 


(5.10) 


E a;,i>[(l Xarc ) EE& x( sin aj k > —, for some k = i,... ,j — l)] 

i=l j> 1 

< p-ZtHee 1 ST' ( 2 fe g7 )° r 7-^ 

few- 1 )' 0 

< CT£ bl ~ s ~ A . 


We can now optimize the parameter <5 setting <5 = 37 2 2 . From equations (5.6) and (5.10) we finally end 
up with 


(5.11) 

which tends to 0 as e —> 0 if 7 G (6/7,1). 


[xrec] — Ce 2 


□ 
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Appendix A. The collision time 


We want to estimate the time spent by a the test particle in the interaction disk associated to the 
central potential of finite range with a uniform magnetic held perpendicular to the plane. Let be e a (f>(r) 
with a £ [0,1/2) the central potential and eB the modulus of the magnetic held. 

The Lagrangian of the system is 

j£?(r, f, 9, 9) = ^ r 2 + ^r 2 0 2 — e a cj)(r) + £-^r 2 9 

We observe that the energy of the system is conserved. Moreover the Lagrangian does not depend on the 
variable 9 , so we obtain the conservation of the conjugate momentum 

—r(r 2 9 + e—r 2 ) = 0. 
dV 2 ' 

Therefore we obtain the following conserved quantities 

r 2 + r 2 9 2 + 2e c V = 2 E, 


B 

r 2 9 + £—r 2 = M, 


and the equations of motion are 


= \1{E-£<*<!>)- 


M 2 e 2 B 2 2 

r 2 4 ' 


eMB 


n d9 M_ -R 

° dt ~ r 2 £ 2 • 


This implies that 


(A.l) 


dt 

dr 

M 

dr 


2 (E - e a cj>(r)) - 


M 2 


e 2 B 2 


-,- 1/2 


■eMB 


M eB 
r 2 2 


I2(E - s a (j)(r)) - -% 2 - ^r 2 + eMB 
We now define the effective potential 

M 2 e 2 B 2 2 eMB 
<t>eff(r) = £ <t>{r) + ^+ g r - 

and we assume that the potential has a short range, i.e. <j>(r) : [0,1] —>• K. and <fi is continuous on [0,1] 
and differentiable on (0,1). 

Take the modulus of the initial velocity to be |u| = 1. When the particle hits the obstacle of radius 
r — 1 the conserved quantities read 

M = p + e 

being p £ [0,1] is the impact parameter. The effective potential is 


<l>eff(r) = £ a (j)(r) 


B 


1 


e 2 B 2 


B\ eB 


P + £ ^) ^2 + — r - ( P + ^ ) — = £ “^( r ) + o 


1 


eB 

~Y 


1 

r- 

r 


By integrating the eciuations of motion we obtain the collision time, namely the time spent inside the 
obstacle: 

2 - 1 - 1/2 


(A.2) 


r = 2 


dr 


eB 


1 - 2 e a <j>{r) - - - — r - - 


1 


where r m i n (the minimum distance from the centre) is the unique zero of the radicand, i.e. 

1 = 2(J) e f f (rmin)i 
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so we can reformulate (A.2) as 


= 72 


dr 


, \/2(0e//( ^ min) 4*e ff ( r )) 

where 2 <j> e ff{r) < 1. The derivative of the effective potential reads 


7 // (r)= £ V(r)-5 3 - 


■ 2 e 2 i? 2 e 2 i? 2 


4r 3 


H- : —r. 


By the mean value theorem we get 


\4*eff (j'min) ~<t>eff{r)\ = \r - r min \\ -<j>' eff {r *)I > |r-r mi „|^ ( inf ^ | -<j>' eff (r)\ 

72 r i 


S {r min ,r) 


and then 


T < 


dr. 


( inf re(r min ,l) I - 7e//7)l) 1/2 Jr min V r ~ 

Since 4>' e ^{r) < 0 for r £ [0,1), then inf re ( rmiTi) i) | — 7 e //( r )I =: K > 0 and it follows easily that 


r < 2 


2(l-r mm )V /2 <2 f2V /2 

K, 


For the corresponding rescaled problem the effective potential reads 


(A.3) 


<P ( eff( r ) = s a (j){r/e) + 


£ P + B _( £ l_ r 


2 r 


with p £ [0,1). In this way one gets —<jr^{r) = ^F(r/e,e) where 
(A-4) F{y, e) = -E a 4>'{y) + (p + B)py~ 3 + ^-( 1 - y 4 ) 

which is positive for y < 1 and uniformly in e. The same argument as before yields the claimed estimate: 


(A.5) t 

where yo = yo{e) is such that 


(2£)V 2 


dy 


{™iye( yo , s )F(y,E)y /2 J yo sjy - y 0 


< Ce 


i=+(- - 7^ + 7 ) ■ 

V 2/o 2 2y 0 / 

We consider now the long range unrescaled potential defined in Assumption B 1), i.e. if e (r) = r_s for 
r < £ 7-1 and ip e (r) = £ _s (t-t) for r > £ 7-1 . The same argument as for the short range case leads to the 
following estimate for the collision time after rescaling: 


(A.6) 

with 


r < 


( 2£) 1 /2 


dy 


< Ce 1 


( inf ye(yo,efo B{y, e )) 1/2 Jyo Vv ~ Vo 
F(y, e) := -ip^y) + E 2h ~ 1] p 2 y~ 3 + B 2 e 3i ~ 2 py~ 3 + B 2 y- 3 £ 47 - 2 ( 1 - y 4 E 4{ - 1 ~^)/A > 0 


for y £ {yojE 1 ) where y 0 = y 0 (e) is such that 1 = 2if e (y 0 ) + 


e 7 1 py 0 1 - 1 y 0 *~ £ yo) 


Appendix B. Cross section 

Proposition B.l. Consider the scattering angle 0{p,E) of a particle with impact parameter p due to a 
uniform magnetic field perpendicular to the plane with modulus eB and due to a radial potential £ a <f>, 
where a > 0 and <f> satisfies assumptions Al, A2, A3. Consider also the scattering angle 9(p,e) associated 
to the same radial potential as before, but without any magnetic field. Then, for e small enough one gets 

(B.l) 0{p,e) = 6{p,e) + O{e). 
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Proof. Following mm we can write the exact formula for both of the scattering angles: 

/■*-»“ {p ' e) du 


(B.2) 


0{p,e) = tt — 2 arcsin p — 2 


where u max {p,£) is the solution of the equation u 


2 

'max 


(B.3) 


0(p, e) = 7r — 2 arcsin p — 2 


r“max(P,E) 


du 


y/l — U 2 — 2 £ a (j)(pu~ 1 ) 
2 = !. while 




where u max (p,£) is the solution of the equation 2e a (j){pu rr f ax ) + u^ ax [ 1 + ^(1 — —)] 2 = 1 . 

Hence, an expansion of 0(p,s) for e small enough yields the claimed asymptotic formula. 


□ 


Proposition B.2. Let 9 be the scattering angle associated to the long range potential 4 >{r) = r s with 
s > 2, 0 £j7 the scattering angle due to a radial potential defined in Assumption B1) and 0e,~j the 
scattering angle due to ip e and to a uniform, constant magnetic field perpendicular to the plane with 
modulus eB. Then one has 

a) 0^7 — > 9 as e —> 0 . 

b) Fe i7 (9) —> r(0) as £ -A 0, where Te,^ (0) is the differential cross section associated to the radial 
potential ip e and the magnetic field, while T( 0 ) is the one associated to the radial potential ’F. 

c) F^((9) < C9~ 1 ~ 1 / s uniformly in £,B. 

Proof, a) Let us now consider the truncated potential >F = r~ s — A~ s with s > 2 for r < A and T = 0 
for r > A with A = e 7_1 and 7 £ (0,1). Take the modulus of the initial velocity of the light particle to 
be |u| = 1 . 

We denote by p the impact parameter (with 0 < p < A) while the scattering angle (that is the angle 
between the ingoing and the outgoing relative velocities) is 


(B-4) 


0e,i(p) = 2 


Pi r/2 


/ arcsin(p/24) 


1 - 


sin/3 


v + sv s 


d/3. 


where v = v{0) such that v 2 + 2{{v/ p) a — A s ) = sin 2 /3 and v = p/r (see Appendix in mr Following 
mu we can write the formula for the scattering angle associated to the potential T and with the 
uniform magnetic field. Due to its invariance under rescaling, the scattering angle associated to the 
equations of motion 2.4 reads 

1 (M-efr 2 )dr 




-2 *eff(r) 


(B.5) 


t/2 


=2 


dfi 


1 - 


(l-^)sin/3 


u + sm s_ 1 M _s + 


e 2 B 2 M 2 
4 u 3 


where 'F e /y(r) = ’F(r) + | (— — | Br) 2 , M is the value of the conserved momentum at the hitting time, 
i.e. M = p + £A 2 ®, and r* is defined as the solution of the equation 2'F e yy(r*) = 1. In the second line 
we made the change r u {3 where u = u(/3, M) = and sin 2 /3 = 2 ^f(M/u). Note that the change 
of variable u —> /3 is well defined because \F(M/u) is non-decreasing when u £ [M/A,M/r*] for e small 
enough. 

From (B.4) and (B.5 1 it is clear that 9 en and 0^ have the same asymptotic behaviour as e approaches 
0. Since 0 £i7 —> 9 , one gets the claim. 

b) The inverse of the differential cross section associated to 9 is 


(B. 6 ) 


y 

2 

I"*/ 2 d/3 sin/3su s 1 

1 + s(s — l)u s 2 p s 

dp 

pS+l J 

arcsin(p/A) ^ P ) 

1 + SV s ~ 2 p~ s 


+ 


2 sp- 2 A 2 ~ s 


1 + 2 sp~ 2 A 2 ~ 
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( R^ 

We want to study the limit of d6g^ /dp. For a mere computational convenience, we prefer to look at 
dO^j /dM which is related to d9ef//dp via 


(B.7) 


From (B.5) one gets 


dd, 


(B) 

e,7 


r m ( B ) dO^ B ^ eBA 2 

w< m >=^ + 4 


sM~ 2 A 2 ~ s + eBA 2 M~ 1 


dM 


A\! i-(f -^) 2 


1 + sM~ 2 A 2 ~ s + § /-B 2 A A M~ 2 




-2 


d/3 sin f3 


— eBM 


aIcsin (^M_eBA ) ( U + SU S ~ 1 M~ S + e2 ^:f P ) 2 

d/3 sin/3 


[sV^AT 8 " 1 - u'( 1 + s(s - 1 )u s - 2 M~ s )] 


f arcsinf — 


M _ eBA \ 
A 2 ) 


u 


s{s + 1)u s ~ 2 M~ s (uM - 1) + 3 u'M 


(B.8) 

where 

(B.9) 


5e 2 B 2 M 3 u' eBM 3eBM 2 3e 2 B 2 M 2 

+ —75 H-—-1 - 


4 u 


2 u 3 


4 m 4 


u = 


su 


s-1 I n sBM \ eBM s + 1 

+ l 1 2u 2 / 2u 


du 1 

dM ~ MS+ 1 ^ l + su s ~ 2 M~ s - ^f /2 


/ D\ 

As for item a), one realizes that ddl tl /dp and d6 en /dp are asymptotically equivalent for any p, thus 
Proposition A.l in m implies that r^(0) —> T(9) for 9 £ (—7r, tt) because its inverse map converges 
everywhere. 

c) From (B .81 for e small enough a tedious expansion gives 
(B.10) 


d9 ( e B) 

> 


C-W (U M =41 ^ 1 

d@e,‘y 

dM 


dp 


dp 


where C > 1 is a constant, 7 Z is bounded in e. The claim follows thanks to Proposition A.l in [DP]. □ 


References 

[BBS] C. Boldrighini, L. A. Bunimovich, and Y. G. Sinai. On the Boltzmann equation for the Lorentz gas. J. Stat. Phys. 
32, 477-501, 1983. 

[BMHH] A. V. Bobylev, F. A. Maa 0 , A. Hansen and E. H. Hauge. Two-Dimensional Magnetotransport according to the 
Classical Lorentz Model. Phys. Rev. Lett., 75, 2, 1995. 

[BMHH1] A. V. Bobylev, F. A. Maa 0 , A. Hansen, and E. H. Hauge. There is more to be learned from the Lorentz model. 
J. Stat. Phys. 87, 1205-1228, 1997. 

[BHPH] A. V. Bobylev, A. Hansen, J. Piasecki and E. H. Hauge. From the Liouville equation to the generalized Boltzmann 
equation for magnetotransport in the 2d Lorentz model. J. Stat. Phys. 102, 1133-1150, 2001. 

[BNP] G. Basile, A. Nota, and M.Pulvirenti. A diffusion limit for a test particle in a random distribution of scatterers. 
J. Stat. Phys., 155, Issue 6, 1087-1111, 2014. 

[BNPP] G. Basile, A. Nota, F. Pezzotti and M.Pulvirenti. Derivation of the Fick’s law for the Lorentz model in a low 
density regime. Commun. Math. Phys., 336, Issue 3, 1607-1636, 2015. 

[DGL] D. Durr, S. Goldstein, and J. Lebowitz. Asymptotic motion of a classical par-ticle in a random potential in two 
dimensions: Landau model. Comm. Math. Phys. 113, 209-230, 1987. 

[DP] L. Desvillettes and M. Pulvirenti. The linear Boltzmann equation for long-range forces: a derivation from particle 
systems. Models Methods Appl. Sci. 9, 1123-1145, 1999. 

[DR] L. Desvillettes and V. Ricci. A rigorous derivation of a linear kinetic equation of Fokker-Planck type in the limit 
of grazing collisions. J. Stat. Phys. 104, 1173—1189, 2001. 

[DR1] L. Desvillettes and V. Ricci. The Boltzmann-grad limit of a stochastic Lorentz Gas in a force field. Bulletin of the 
Institute of Mathematics Academia Sinica (New Series) 2, 2, 637-648, 2007. 

[DR2] L. Desvillettes and V. Ricci. Non-Markovianity of the Boltzmann-Grad limit of a system of random obstacles in a 
given force field. Bull. Sci. Math. 128, 1, 39-46, 2004. 

[G] G. Gallavotti. Grad-Boltzmann limit and Lorentz’s Gas. Statistical Mechanics. A short treatise. Appendix 1.A2. 
Springer, Berlin, 1999. 

[H] E. H. Hauge. What can one learn from the Lorentz model? Transport Phenomena, G. Kircenow and J. Marro, 
eds., Lecture Notes in Physics, Vol.12, Springer-Verlag, Berlin, 1974. 

[K] K. Kirkpatrick. Rigorous derivation of the Landau Equation in the weak coupling limit. Comm. Pure Appl. Math. 
8, 1895-1916, 2009. 




































22 


M. MARCOZZI AND A. NOTA 


[KP] H. Kesten and G. Papanicolaou. A limit theorem for stochastic acceleration. Commun. Math. Phys. 78, 19-63, 
1981. 

[KS] A. Kuzmany and H. Spohn. Magnetotransport in the two-dimensional Lorentz gas. Phys. Rev. E , 57, 5, 1998 . 

[L] Landau, L.D., Lifshitz, E.M.: Mechanics. Course of Theoretical Physics, vol.l. Pergamon press, Oxford, 1960. 

[S] H. Spohn. The Lorentz flight process converges to a random flight process. Comm. Math. Phys. 60, 277-290, 1978. 

(Matteo Marcozzi) University of Helsinki, Department of Mathematics and Statistics FI-00014 Helsingin 
yliopisto, Finland 

E-mail address , M. Marcozzi: matteo.marcozzi@helsinki.fi 

(Alessia Nota) Institute for Applied Mathematics, University of Bonn, Endenicher Allee 60, D-53115 Bonn, 
Germany 

E-mail address, A. Nota: nota@iam.uni-bonn.de 



